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Abstract 

We introduce and study a deterministic lattice model describing the motion of an infinite 
system of oppositely charged particles under the action of a constant electric field. As an 
application this model represents a traffic ffow of cars moving in opposite directions along a 
narrow road. Our main results concern the Fundamental diagram of the system describing the 
dependence of average particle velocities on their densities and the Phase diagram describing 
the partition of the space of particle configurations into regions having different qualitative 
properties, which we identify with free, jammed and hysteresis phases. 
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1 Introduction 

It is well known that far from equilibrium statistical physics systems may demonstrate a very 
complex behavior even on the level of their stationary states. In general, such non-equihbrium 
stationary (steady) states depend sensitively on the details of the microscopic dynamics which 
leads to a very diverse organization at the macroscopic level. An important class of such systems 
is represented by driven diffusive systems (DDS), describing the behavior of interacting diffusing 
particles driven into selected spatial directions by an external force. Models of this sort are 
widely used to describe vehicular and pedestrian traffic [2T], water droplets in microemulsions 
with distinct charges [2j, and numerous biological problems from molecular motors to protein 
synthesis (see e.g. [23j). 

Basic ideas used in these models were introduced originally in [19j as a simple modification 
of the Ising lattice gas. Prom that time tens of publications have been devoted to this subject. 
From the point of view of mathematics these processes are close relatives to asymmetric exclusion 
processes and topological Markov chains. 

The simplest model of this sort can be described as follows. Consider a finite one-dimensional 
lattice with periodic boundary conditions (i.e. a ring). Each site of the lattice is either vacant 
(•) or occupied by a positive (+) or negative (-) particle. Positive particles are driven to the 
right and negative particles are driven to the left. The dynamics is random-sequential, i.e. at 
each time step a pair of nearest neighbor sites is selected at random and particles or vacancies 
located at these sites are exchanged according to the following set of probabilities: 

Vi+- ^ +-)=/?, Vi-+ ^ +-) = 0, Vi+ ■ ^ ■+)=Vi > -■) = !. 
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Here the parameter < /3 < 1 plays the role the time delay required to exchange neighboring 
positive and negative particles H — . Numerical studies and mean-field type approximations of 
this model and its two-lane generalizations have been a subject of a large number of publications 
(see e.g. [Il[l3l[9l[IS[I7lligiIIl[l2l|20l[ini|2l[l5]and further references therein) . 

Our aim is to introduce a deterministic version of a DDS and to study rigorously its limit 
dynamics (as time goes to infinity) on the infinite integer lattice Z. In order to do this we first 
pass from the random-sequential to synchronous movement of particles (when all particles are 
trying to move simultaneously), and then construct deterministic exchange rules for particles of 
different signs. The simplest way to achieve the second goal is to replace the random exchange 
with the introduction of the waiting time t = 1/ (3 before the exchange. (Indeed, the average 
waiting time in the random exchange version is exactly 1//3.) Further on we assume that the 
parameter r takes only positive integer values. 

To define the synchronous movement of particles one needs to resolve the problem with 
the presence of triples -|- • — (positive and negative particles separated by a single vacancy) 
additionally to couples H — considered in the asynchronous case. In the case of a triple -|- • — 
both positive and negative particles are trying simultaneously to exchange their positions with 
the same vacancy (which cannot happen under the asynchronous updating). We assume that 
in this case the waiting time before the particles exchange their positions is equal to 
The situations H — and -|- • — we call short and long interactions^ while the situations -|— |- and 

when a particle is blocked by another particle of the same sign we call simple interactions. 

Using the notion of interactions the dynamics can be described as follows. All particles in a 
configuration are trying to move in the direction corresponding to its sign by one position. If 
the particle does not interact with others then it makes a move. If the short/long interaction 
takes place then the particle waits for r or r -|- 1 time steps (depending on the interaction type) 
before to start moving, while in the case of the simple interaction the particle simply preserves 
its position. As an example consider the dynamics of a spatially periodic configuration with 
r = 2 (only the main period of length 17 containing 4 positive (+) and 2 negative (~) particles 
and 11 vacancies ( . ) is shown): 
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Here the first two positive particles and the first negative particle move freely until the 
short interactions with the 1st negative particle consecutively take place at time t=2 and t=4. 
Then after 2 time steps the positive and negative particles exchange their positions. The 3d 
positive particle at the 1st time step makes a simple interaction with the 4th positive particle. 
The 4th positive particle starting from t=l is making a long interaction with the 2nd negative 



^ One might argue that after the first time step both particles move to the vacant site and then use the 
standard waiting time r to make further exchanges. 



Figure 1: (left) Fundamental diagram: dependence of the limit average velocity V on the density 
p of positive particles with the fixed density p < l/(r + 2) of negative particles. The average 
velocity V of negative particles is indicated by a thin line. Values PcP'c indicate boundaries of 
the hysteresis phase. 

(right) Phase diagram: A C Free+ n Free!^, B C Free+ , C C Free^f . The region H := {0 < 
p + p < 1}\{AL) B L) C) between thick and thin lines belongs to the hysteresis phase. 

particle during three time steps, then exchange their positions, etc. Note that due to the spatial 
periodicity after the 5th time step a positive particle "moves" from the 17th to the 1st position, 
while after the 12th time step a negative particle "moves" from the 1st to the 17th position. 

In fact, this model is not Markov since the information about present positions of particles 
alone does not allow to predict the future motion of particles (due to the time delays during 
particle interactions) and to make the model precise one needs to introduce the notion of particles 
states which will be done in Section [2j 

The analogy between the classical gas-liquid transition and stationary states in interacting 
particle systems is well known and discussed at length in the literature. The gaseous or laminar 
phase corresponds to the situation when particles almost never interact with each other and thus 
move with constant velocities. Typically this happens at low particle densities. In the liquid 
phase (which occurs normally at high densities) the particles are located close to each other and 
mutual interactions organize them into clusters (traffic jams). Additionally there might be an 
intermediate or hysteresis phase when both types of particle behavior coexist. 

In our case the classical laminar phase cannot take place since particles moving in opposite 
directions interact inevitably. To overcome this difficulty we say that a particle is becoming free 
starting from a certain time to ^ if for Vt > to it does not take part in simple interactions. 
In other words, after an initial transient period of length to the particle stops interacting with 
particles of the same sign. Then instead of the laminar phase we consider a pair of eventually 
free positive and negative phases defined as subsets Free^ of the set of all particle configurations 
as follows. A configuration belongs to FreeJ^ if each its positive particle becomes free after a 
finite transient period, which might depend on the particle. The set Free^^ is defined in the 
same way but for negative particles. Similarly instead of the liquid phase we consider two 
eventually jammed (notation Jam^) phases having the property that clusters of particles of the 
corresponding sign are always present after some transient period. The cases when the length of 
the transient period is equal to zero we denote by Free-t and Jam-t respectively. Note that the 
sets of eventually free and eventually jammed configurations of different signs have nonempty 
intersections and that it is possible that an initial configuration has no interacting particles but 
under dynamics they will start interacting and form "jams". 

Our main results are formulated in the following two theorems. We start with qualitative 
results justifying our discussion of phase transitions and giving an inner characterization of 
individual configurations belonging to each phase. To this end in Section [3] we introduce the 



notion of signed (positive and negative) (proto)clusterqj of particles. By the life-time of a 
(proto) cluster one means the duration of time until it ceases to exist. Exact definitions and 
discussion see in Sections I2|3I 



Theorem 1.1 Let for a configuration x the densities of positive and negative particles {p{x), p{x)) 
be well defined. If 

(a) (r + 2)p{x) < 1 + (r - l)pix) then x G Free!^ 

(b) (r + 2)p{x) < 1 + (r - l)p{x) then x G Pree^ 

(c) (r + l)p{x) > 1 + (r - l)p{x) then x £ Jam° 

(d) (r + l)p{x) > 1 + (r — l)p{x) then x G Jam^r 
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To formulate quantitative results consider a partition of the triangle A := {{p,p) 
p + p < I, p, p > 0} (describing all possible pairs of densities) made by 4 straight lines: 
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{r + l)p 
(see Fig [1] (right)). Denote 
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(r + 2)p < 1 + (r - l)p, (r + 2)p < 1 + (r 
(r + 2)p< 1 + (r - l)p, (r + l)p > 1 + (r 
(t + 1)p> 1 + (t-1)p, (r + 2)p<l + (r 
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Theorem O implies that A C Free+ n Free^f , B C Free+ n Jam°^, C C Free^f n Jam5 



< 



Theorem 1.2 Let the densities of both positive and negative particles p,p in the initial config- 
uration be well defined. Then the corresponding average particle velocities V, V are well defined 
as well and depend only on the particle densities according to the following relations: 

np,P) = lZ-lKJr np,P) = \t[:zi\tJ) ^f{p,p)eA 

V{p,p) = ^, y(p,p) = i (i-i) if{p,p)eB 
Vip,p) = ^ Vip,p) = ^ ifip,p)eC 

Note that Theorem 11.11 gives a qualitative description of the asymptotic dynamics, while The- 
orem 11.21 describes it in quantitative terms. It is worth mention also that without a systematic 
preliminary numerical modelling which allow us to understand qualitatively the local structure 
of the particle flow even the formulation of results proven in this paper would be impossible. 

Proofs of these results ideologically are based on the machinery of the analysis of life-times 
of particle clusters developed in ^ [SJ [6] where deterministic interacting particle systems with a 
single type of particles were studied. The presence of particles moving in the opposite direction 
together with two types of interactions (short and long) complicates significantly the behavior 
of the system. In particular, new clusters may be born in a free flow of particles, and the 
asymptotic dynamics depends on two parameters (densities of positive and negative particles). 
Therefore the technics has been changed a lot and yet we are able to give only lower and upper 
estimates of the life-times. Nevertheless these estimates allow to find exact boundaries of all 
phases present in the model. 

^Roughly speaking a protocluster is a collection of particles of the same sign which will form a true cluster 
in future. A configuration x £ Jam^ may contain no infinite life-time clusters but in that case there are infinite 
life-time protoclusters. 



The paper is organized as follows. In Section [2] we give the formal description of the model 
and main statistical quantities under study: particle densities, average velocities, etc. In the 
absence of clusters the dynamics of particles is trivial (except from interactions between par- 
ticles of opposite signs). Therefore to study the model we need to analyze the dynamics of 
clusters of particles and their "life-times" . Exact definitions of these objects and corresponding 
mathematical results are discussed in Section [3l Duality relations between positive and negative 
particles allow us to calculate in Section [4] average particle velocities under the assumption that 
they are well defined. The latter is connected to the proof of our main results formulated in 
the Introduction and given in Section [5l Section [6] is dedicated to the analysis of the region of 
the Phase diagram where different phases coexist. In Section [7] we apply our results to study a 
model of an active tracer moving with or against the particle flow, and Section [8] is dedicated to 
the generalization of qualitative results for the case of configurations for which particle densities 
are not well defined. 

2 The model 

The main disadvantage of the model discussed in the Introduction is that the information about 
the present positions of particles alone without the knowledge for how long currently occurring 
interactions already take place does not allow to define the future motion of particles. In order to 
overcome this difficulty we introduce the notion of a state of a particle which takes into account 
the complete information about the occurring interactions. Let us give formal definitions. 

By a configuration we mean a bi-infinite sequence x = {xi}^^ with elements from the 
alphabet {— r — 1, — r, . . . , —1, 0, 1, . . . , r, r-|- 1}. Positive entries correspond to positive particles, 
negative entries to negative ones, while zero entries correspond to vacancies. Non zero entries 
will be referred as states of particles located at corresponding sites. The states will be used to 
take into account the delays during interactions of positive and negative particles. Thus the 
largest value of the state t + 1 corresponds to the long interaction. 

The set of admissible configurations X consists only of configurations x satisfying the con- 
dition that for each i E Z 

• if > 1 then either Xj+i = — Xj or Xj+i = and Xj_|_2 = —Xi, additionally Xi = t + 1 
implies Xj+i = 0; 

• if < —1 then either Xi-i = —Xi or Xi-i = and Xj_2 = — Xj, additionally Xj = — r — 1 
implies Xj_i = 0. 

These conditions imply restrictions only to positions of interacting particles having states greater 
than one on modulus. For example, the configuration . . . 110101100020200 ... is admissible (here 
1 and 2 stand for —1 and —2), while . . . 110101100020100 ... is not, since the last two particles 
separated by a single vacancy 201 are supposed to be mutually interacting (according to their 
positions) but their states differ on modulus. 

By X* := r*x, t G U {0} we denote the state of the configuration x at time t, assuming 
that the initial state x'' is given by x. Here T : X ^ X is the map describing the dynamics 
which we define on the level of individual particles in the configuration x £ X in the following 
three steps: 

1. First consider sites i,i' with |i' — z| < 2 containing mutually interacting particles with 
X* > 1 and X*, = — x* and set xf^^ := x* -|- 1, x*,"^^ := x*, — 1. Then if — i = 1 
(short interaction) and x*^^ > r or if i' — i = 2 (long interaction) and x*"*"^ > t + 1 set 
xr^:=-l, x*+^:=l. 



2. Then consider the sites i with x\ = 1. 

(a) if x-_^i = and x\j^2 / "1 then set := 0, x*+} := 1; 

(b) if = and x*+2 = -1 then set xf-^ := 2, x-+2 := -2; 

(c) otherwise if x\j^^ = —1 then set x*"*"^ := 2, := —2. 

3. It remains to consider the sites i with x\ = —1 which were not taken into account during 
the step 2. If x\_i = 1 set := —1, x*"*"^ := 0; otherwise do nothing. 

In words, if a particle is not interacting it simply moves by one position in the direction 
corresponding to its sign (rule 2. a and 3). In case of the simple interaction (the particle is 
blocked by another particle of the same sign) it does not move an does not change its state. If 
the short/long interaction takes place the particle preserves its position but its state changes 
by ±1 depending on the particle sign (see rule 1) until it reaches on modulus the value r (in 
the case of the short interaction) or r + 1 (in the case of the long interaction). After that the 
particles get the states ±1 (preserving original signs) and exchange their positions. The rules 
2.b and 2.c take care about the initial stage of interactions. 

Using the notation 1 = —1,2 = —2,3 = —3 we can rewrite the example of dynamics of a 
spatially periodic configuration (see the formal definition below) with r = 2 described in the 
Introduction as follows: 
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Invariance of the set of admissible configurations follows immediately from the definition 
of the map T. Observe also that the restriction of the map T to the subset of admissible 
configurations containing only entries from the alphabet {0, 1} coincides with the classical Nagel- 
Schreckenberg traffic flow model (see e.g. |221 [5]). 

For a configuration x € X by the density of positive particles p{x, I) in a finite lattice 
segmenlH / = [n, m] := {i € Z : n < i < m} we mean the number of positive particles from 
the configuration x located in / divided by the total number of sites in / (notation |/|), and by 



^We shall use also the notation x[n, m] for lattice segments to specify the configuration x. 



p{x, I) the corresponding value for the negative particles. If for any sequence of nested finite 
lattice segments {In} with oo the limits 

p{x) := lim p{x,In), p{x) := lim p{xjn) 

are well defined and do not depend on {In} we call p{x) the density of positive particles and 
p(x) the density of negative particles in the configuration x € X. Otherwise one considers upper 
and lower particle densities p^{x), p^{x). 

Lemma 2.1 Particle densities are conserved under dynamics, i.e. p^{x^), p^{x^) do not depend 
on t < oo. 

Proof. For a given lattice segment / G Z the number of particles from the configuration 
X* G X which can leave it during the next time step cannot exceed 2 and the number of particles 
which can enter this segment also cannot exceed 2. A close look shows that the total number of 
particles that can leave or enter this segment cannot exceed 2, because if a particle leaves the 
segment through one of its ends no other particle can enter through the same end. Therefore 

\p{x\l)-p{x'+\l)\-\I\<2 

which implies the claim in the case of positive particles. The proof in the case of negative 
particles is exactly the same. □ 

A configuration x £ X is said to be spatially periodic with period L € Z if = Xj+i Vz € Z. 
Such configurations represent an important class of admissible configurations for which both 
positive and negative particle densities are well defined. A direct check demonstrates that the 
spatial periodicity and its periocl^ are preserved under dynamics. Therefore the dynamics of 
spatially periodic configurations of period L is equivalent to the dynamics on the finite segment 
of the integer lattice of length L with periodic boundary conditions (i.e. to the only situation 
suitable for numerical modelling). The total number of admissible configurations on this lattice 
segment cannot exceed L^'^^^, which implies eventual periodicity in timqj of the dynamics. Thus 
a trajectory {r*x}t>o of the original infinite system with a spatially periodic initial configuration 
X € X is eventually periodic in time as well. 

We say that particles of the same sign located at sites i' < i are consecutive if all sites 
between them are either vacant or occupied by particles of the opposite sign. Consider a pair 
of consecutive positive particles. The following result shows that the distance between these 
particles (calculated as i — i' — 1) can shrink to zero with time but cannot be enlarged much. 
(Compare with the incompressibility property of a conventional fiuid!) 

Lemma 2.2 Let X be a collection of admissible configurations x having a pair of consecutive 
positive particles at sites i' < i. Denote by f'(x*),i(x*) positions of these particles at time t in 
the configuration x* and by D(x*) := i(x*) — i'(x*) — 1 the distance between them. Then 

= inf liminfL)(x*) < sup lim sup D(x*) < 2r(z -i'). (2.1) 

Proof. The lower estimate follows from the observation that, independently on the initial 
distance between the particles, there might be a large enough "jam" of positive particles ahead 
of them which will stop the leading particle and allow the rear one to catch up with it. 



*The minimal spatial period may decrease under dynamics. 

trajectory becomes periodic in time after some finite transient period. 



To prove the upper estimate observe that without interactions with negative particles the 
distance D{x^) cannot grow in time. Note also that a finite number of negative particles located 
initially in the segment [i' , i] can give only a constant contribution to the variation of the distance 
between the particles. Therefore it is enough to consider only the case when initially there are 
only vacancies in the segment [i' + l,i — 1]. Denote by ti the first moment of time when the 
leading particle meets a negative one, by t[ the duration of the interaction between them (which 
might take values r or t + 1), by t2 the time between the end of this interaction and the moment 
when the rear particle meets with the same negative particle, and by ^2 ~ the duration of the 
latter interaction. Then we have 

= max{f'(x*i) +t;,i(x*i) - 1}, 
=i(x*i) + l + (t;-r); 

< + 1 + {t[-T)+t2; 

i'^^h+t[+t,+t',-^ ^ .;^^t,+t[+t,^ + 1 + (4 - r), 

i(x*l+*'l+*^+*^) < + 1 + -r)+t2+t'2. 

Therefore 

L)(x*i+*'i) = i(x*i+*i) - i'(x*i+*'i) - 1 

= max{i(x*i) - - r + 1 - 1,2} < max{Z)(x*i), 2}, 

< i(x*i) + (4 - r) + t2 - max{i'(x*i) + t[,i{x*') - 1} - t2 

< max{i(x*i) - i{x^') -t -2,2} < max{Z)(x*i), 2}, 

£)(a;*i+*'l+*2+4) = ^(-^il+i'i+*2+t2) _ i'(x*l+*i+*2+*'2) - 1 

< i(x*i) + (t'l -T)+t2 + t'2 - max{i'(x*i) + t[,i{x^') - 1} - t2 - I - {t^ - t) 
= i(x*i) + t[-l- max{i'(x*i) + t[,i{x^^) - 2} 

< max{i(x*i) - i'(x*i) - 1,*^} < max{D(x*i), r + 1}. 

This implies the result since the contribution from a single negative particle located initially in 
the segment cannot exceed r + 1 and the number of such particles cannot be larger than 
i-i' -1. □ 



Remark 2.3 Result of Lemma 12.21 is not obvious and is not what one expects here. Indeed, 
due to the presence of two types of interactions: short and long, one would expect that if the 
leading particle in the pair takes part only in short interactions while another one takes part 
only in long interactions, then the distance between them grows linearly in time. Lemma 12.21 
demonstrates that this is not the case. An important property that we use here is that the 
difference between the duration of long and short interactions is exactly equal to one. If this 
would not be the case the distance between the particles indeed might grow with time. 

Let us introduce now the notion of the average velocity of a particle. For a configuration 
X denote by the location of a certain particle at time t > located originally at site i = i^. 
By a finite time velocity of this particle we mean V{x,i,t) := J |i* If the limit V{x,i) : = 

lim^^^oo Vix, i, t) is well defined we call it the average velocity of the particle. To distinguish 
between the positive and negative particles in the latter case we use the notation V{x,i). 

Our aim now is to show that the average velocity does not depend on the choice of a particle. 

Lemma 2.4 Let x he an admissible configuration and assume that for a positive particle origi- 
nally located at site i the average velocity V{x,i) is well defined. Then for any positive particle 
in the configuration x the average velocity is well defined and coincides with V{x,i). 



Proof. Denote by i the location at time t = of the next positive particle located to the right 
from i. As usual we denote by t > positions of these particles at time t > 0. Then we 

have 

Applying Lemma 12.21 according to which 1 < i* — i* < 2r|i — i|, we get 

\V{x,i,t)-V{x,i,t)\ < J + ^ |i° < \i-i\ *-^0. 

Thus V{x, i) = V{x, i). Using the same argument one extends this result to neighboring positive 
particles, and repeating it to all positive particles in the configuration. □ 



Corollary 2.5 The invariance of velocities of negative particles can be proven along the same 
argument. Using these results we may drop the dependence on the index i in the definition of 
the average velocity. 

3 Dynamics of (proto)clusters and their basins of attraction 
(BA) 

For a collection of consecutive particles of the same sign ordered with respect to their sign (i.e. 
according to their positions in the case of the positive sign, and in the opposite way in the case 
of the negative sign) we call rear/leading the first/last particle with respect to this order. By a 
cluster (of particles) we mean a segment J = x[m' , m], m' < m of the configuration x consisting 
only of particles of the same sign (called the body of the cluster) and in which the leading particle 
is interacting with a particle of the opposite sign (called the root of the cluster) and the site 
next to the rear particle is vacant or is occupied by a particle of the same sign as the root. The 
sign of a cluster is identified with the sign of particles in its body. Depending on the type of 
the interaction between the leading particle and the root we say that the cluster has a short or 
long root. The right hand side of the diagram below shows positive clusters with short and long 
roots, while the left hand side demonstrates how these clusters may be created: 

'■■■++, ■■■ 

body short root 
■■■++,■ ••• 

body long root 

In distinction to other deterministic transport flow models (see e.g. OHlEllG]) in the model 
under consideration new clusters might be created even in the absence of clusters in the initial 
configuration. 

In order to control events of this sort we introduce the notion of a protocluste^ as a pair of 
consecutive particles of the same sign and a particle of the opposite sign satisfying the property 
that these three particles will form a true cluster and the moment when the leading of these two 
particles will start interacting with the rear one is the moment of the creation of a new cluster. 
In other words, they do not join an existing cluster. The minimal lattice segment containing 
this pair of particles is called the body of a protocluster and the particle of the opposite sign 
with which they will form the cluster - its root. The sign of a protocluster is again identified 



In the analysis of the free phase we will need to prove that the conditions of Lemma 13. 31 imply that a flow of 
noninteracting at time t > positive particles remain free under dynamics. 



with the sign of particles in its body. Example of a positive protocluster and clusters which 
it may form are shown in the diagram above. Note that between the body and the root of a 
protocluster there might be some other particles of the same sign as the body, which even might 
form a cluster with the root but will stop interacting with the root particle before the leading 
particle of the protocluster will reach them. 

Remark 3.1 A simple calculation shows that particles of the same sign as the root might be 
present between the body of the protocluster and its root only if they will make short interactions 
with the particles from the body while with the root they will make long interactions. Indeed, 
only in this case a new cluster can be created. 

The boundaries of the body of a (proto)cluster might change with time when the particles are 
moving (in the case of a protocluster) or the leading particle leaves it while some particles join 
the cluster from the other side. At time t > by J* we denote the segment of the configuration 
corresponding to the body of the (proto)cluster originally located at J := J*^. If t = we 
drop the time index to simplify the notation. The duration of time during which the cluster 
exists (i.e. the number of particles in its body is not smaller than two and the leading one is 
interacting with the same root particle) we call the life-time of the cluster J* and denote by 
LT{J). By the life-time of a protocluster we mean the duration of time until the cluster formed 
by particles from the body of the protocluster ceases to exist. If during its life-time the body 
of a (proto) cluster consists of only two particles (being present through its entire life-time) we 
call it trivial and nontrivial in the opposite case. Clearly the life-time of a trivial (proto) cluster 
is finite and depends only on the distance to the root. Therefore we shall be interested only 
in the case of nontrivial (proto) clusters. Observe that the situation described in Remark 13.11 
corresponds to the trivial protocluster. 

'Attracting' preceding particles, a (proto)cluster J plays a role similar to an attractor in 
dynamical systems theory. Therefore it is reasonable to study it in a similar way and to in- 
troduce the notion of its basin of attraction (notation BA(J)), by which we mean the minimal 
segment of the configuration x containing all sites from where particles may eventually join the 
(proto)cluster during its life-time. 

Due to the duality between positive and negative particles it is enough to consider only 
positive (proto)clusters. Let J = x[m',m], m' < m he the body of a positive (proto) cluster 
having its root at site m" > m. To take into account the dynamics we use the notation 
J* = x^[m' {t),m{t)] and m"{t) > m{t) for the corresponding objects at time < t <LT(J). We 
introduce two functionals depending on an integer parameter k < m': 

m m 

W{x[k,m]) :=q{xm) + ^w{xi), W'{x[k,m]) := q{xm) + ^w'{xi), (3.1) 

i=k i=k 

where 

{T + iifz>o (t a z>0 

—T if z < , w'{z) := < — r if z < , 

—1 otherwise [^ — 1 otherwise 

and q{z) is defined as the amount of time until the leading particle leaves the cluster, or q{z) := r 
in the case of a true protocluster. Thus q{z) = r — z + 1 for the cluster with the short root, and 
g(z) = r — z + 2 for the cluster with the long root. 

Denote by j{t) the largest integer non-exceeding m'{t) for which W{x*[j{t),m{t)]) = 
or W{x^\j{t) — l,m{t)]) < 0, and by j'{t) the largest integer non-exceeding m'{t) for which 
W{x^[j'{t),m{t)]) = or W{x^[j'{t) - l,m{t)]) < 0. Note that both j{t) and j'{t) might take 
an infinite value. 

The following result allows to control the left boundary of a BA. 



Theorem 3.1 x[j' {t),m{t)] CBA{J^) C x[j{t),m{t)] andj{t+l) > j{t) + l, j'{t + l) = j'{t) + l. 



The proof of this Theorem wih be given through a series of Lemmas, but before to do this 
consider an example demonstrating that the boundaries j, j' might differ a lot and the quantity 
j' — j even might take an infinite value. Consider a spatially periodic configuration x with the 
main period containing three particles and three vacancies: + • • + • — . These three particles 
represent a protocluster for which j(0) = —oo, j'(0) = i — 2 where i stands for the position of 
the 1st positive particle. 

Lemma 3.2 The life-time of a cluster of positive particles is equal to the number of positive 
particles in its BA multiplied by r minus the state of the leading particle plus one in the case of 
the short root or plus two in the case of the long root. 

Proof. The term "the number of positive particles in its BA multiplied by r" describes the 
amount of time before all particles in the BA will leave the cluster, assuming that the initial 
state of the leading particle is one. Indeed, during the life-time of a cluster the current leading 
particle is making continuously the short interaction with the root particle. After this the role 
of the leading particle is going to the next particle in the original BA. The term "the state of 
the leading particle plus one ..." takes into account the actual state of the leading particle and 
the type of interactions with the root. □ 



Lemma 3.3 BA(J)C. x[j,rn\ and the life-time of the cluster J can be estimated from above as 
the number of positive particles in x\j,m] multiplied by (r + 1). 

Proof. Consider the segment / := x[j,m], set n := m — j + 1, and denote by N the number 
of positive particles in this segment (which cannot be smaller than 2 by the definition of the 
cluster) and by K the number of vacancies in it. Then the number of negative particles in this 
segment is equal ton+1 — A^ — i^. Assume contrary to our claim that there exists a particle 
^ located initially to the left from the site j which joins the cluster J := x[m',m] during its 
life-time. Since both vacancies and negative particles are moving to the left (opposite to the 
movement of positive particles) the particle ^ must interact with all these vacancies and negative 
particles before to join the cluster J. The amount of time necessary in order to do this is at least 
t' := K -\- (n -\- 1 — N — K)t. By the definition of the functional W{x[j, m]) and the parameter 
j we have 

N{t + l)-K-{n + l-N- K)t < 0. 

Using this inequality we estimate the time t' from above as N{t + 1), which is larger than the 
upper estimate Nt -|- 1 of the amount time during which all positive particles will join the 
cluster J and then leave it. Therefore the particle ^ cannot join the cluster J in time. We came 
to the contradiction. □ 



Lemma 3.4 Let J be a nontrivial protocluster, then BA(J^) C x[j{t),m{t)]. 

Proof. Denote by ^(i) the position at time t > of a particle being originally at time t = 
to the left from the site j(0). Our aim is to show that ^{t) < j{t) for any moment of time t > 
until the cluster is formed. Whence it is formed we can use the result of Lemma 13.31 

Observe that during the initial part of the life-time of a protocluster the leading particle 
moves at rate one exchanging positions with vacancies until at time ti it starts interacting with 
the root. After this whence another particle from the body of the protocluster reaches the 



leading one at time t2 > ti the cluster is formed and one can apply Lemma [3.3l to get the result. 
Therefore we are concerned only with these two initial parts of the dynamics of a protocluster 
mentioned above. 

During the time < t < ti the leading particle is moving freely exchanging its position with 
vacancies. Therefore 



because a new vacancy is taken into account while all already present particles and vacancies 
contribute the same values as before. 

Observe that the site j{t) cannot be occupied by a positive particle, otherwise 



which contradicts to the definition of j{t). Thus the boundary j{t) moves at least at rate one. 
On the other hand, the particle ^ also cannot move faster than at rate one. Therefore ^(0) < j(0) 
implies that ^(t) < j{t) during the initial part of the life-time of a protocluster. 

During the time ti < t < t2 the leading particle stays at the same site and thus there are no 
new vacancies. However, the state of the leading particle is increasing at rate one which plays 
exactly the same role and leads to the same result as in the previous case. □ 

To deal with the jammed phase we need to consider the lower bound for the BA. To this end 
we apply the functional W instead of W to show that the corresponding boundary is located 
inside of the BA and that all particles from the corresponding segment must join the cluster 
during its life-time. 

Lemma 3.5 5^(J*) D x[j'it),m{t)]. 

Proof. The argument here is very similar to the one used in the proof of Lemma 13.41 with 
the only difference that we need to show now that for any particle ^{t) initially located in the 
segment x[j(0), m(0)] we have ^(t) > j{t). In other words, that the boundary cannot outran 
the particle ^{t). 

The state of the current leading particle of the cluster is increasing by one until it will 
exchange positions with the root particle. Therefore the value of the functional W' is decreasing 
at rate one. Observe now that the site j'{t) cannot be occupied by a positive particle (otherwise 
W{x^[j' {t),m{t)]) > t) and hence W'{x^[j'{t) — l,m(t)]) > which contradicts to the definition 
of j'(t)). Thus the boundary j'{t) moves at least by rate one. On the other hand, the particle 
^ also cannot move faster than at rate one. Therefore j{t) being initially smaller than ^(t) may 
outran the latter only if at some time t' it makes a jump longer than one. This can happen only 
if there is i such that j(t' — 1) < i < m{t' — 1) and m(t' — 1)]) = 1. In this case one 

might expect that on the next time step we have W'{x^ [i,m{t')]) = 0. However this happens 
only if the segment x* [i — 1, i] at time t' — 1 is occupied only by vacancies, which in turn means 
that W'{x^ [i — l,m{t' — 1)]) = and hence i — 1 = j{t'). Thus long boundary jumps cannot 
take place. □ 

This finishes the proof of Theorem 13. 1[ 

The functional W' is everywhere positive in the case of an infinite BA and is used to calculate 
the second critical density corresponding to the beginning of the jammed phase: 



W{x'+\j{t)Mt + m = W{x'[j{t),m{t)]) - 1 



W{x^[j{t),m{t)]) > r + 1 



and hence 



W{x^[j{t) - l,m{t)]) > 



pr — pr — {1 — p — p) > 




l + p(r-l) 



Pc = 



1 -I- T 



In the "intermediate phase" (see Section [6]) and during the transient period of the free phase 
there are no infinite hfe-time clusters but instead there might be recurrent ones, which are 
appearing and disappearing again (in the intermediate phase) and finite hfe-time clusters even 
in the eventually free phase. To calculate the corresponding critical density pc we need take 
into account that the duration of the long interaction is r + 1 instead of r. Therefore it can be 
calculated using the functional W in which the weight of a positive particle is (r + 1) but the 
weight of a negative particle is — r. The condition l^(x[m — n, m]) < for some arbitrary large 
n gives 

p{r + l)-pT-{l-p-p)>0 =^ p(2 + t)-p{t-1)>1 =^ = i±|il^. 

4 Duality relations and calculation of average velocities 

In this Section we discuss connections between dynamics of particles of different types. 

A configuration (— x) € X is said to be dual to x. This is indeed the case since the state of 
each particle is inverted from positive to negative and vice versa. It is straightforward to check 
that the dynamics of configurations x* and — x* are exactly the same except for the direction 
of the particles movement. This type of relations we call duality and they allow to simplify 
significantly the analysis of such systems. In particular we shall mainly discuss properties of 
positive particles having in mind that the corresponding properties of negative particles can be 
obtained by duality relations. 

For a configuration x & X with densities of positive and negative particles p{x) and p{x) 
respectively denote by V{p, p) the average velocity (if it exists) of positive particles and by V{p, p) 
the modulus of the average velocity of negative particles. By the duality between positive and 
negative particles we have V{p,p) = V{p,p). 

It turns out that the assumption about the existence of average velocities together with some 
information about the structure of configurations allows to calculate the average velocities in 
terms of particle densities. 

Lemma 4.1 Assume that in a configuration x the densities and average velocities of both types 
of particles are well defined and that positive particles move freely (i.e. except from the interac- 
tions with negative particles). Then 

V{p, p) = 1 - (r - l)p (Vip, p) + V{p, p)) (4.1) 

Proof. During the time t > a positive particle will meet {y{p, p) + V{p, p)) •tp + o{t) negative 
ones. The notation o{t) stands for a term which grows slower than t, i.e. o{t)/t g. Thus 
taking into account that according to the assumption positive particles never wait in clusters 
(which otherwise will give an additional contribution) we have: 

V{p, p).t = t-{T- l){V{p, p) + V{p, p)) ■ tp + o{t). 

Dividing both sides of this relation by t and passing to the limit as t — > cxd we get the result. □ 



Lemma 4.2 Assume additionally to the assumptions of Lemma \4.1\ that negative particles are 
free also. Then 

,,, ^> 1 + (t - l){p - p) 



Proof. By Lemma l4.ll we have 

v{prp) = i-{T-i)p{v{prp) + v{prp)) 
v{prp) =i-{T-i)p{v{p,p) + v{p,~p)). 

Denote v = V{p, p) + V{p, p). Then adding the above relations we get v = 2 — {t — l){p + p)v, 
which yields 



1 + (t-1)(p + p) 

and hence 

l + ir-l)ip + p) 

In particular, 

V{p,p) ^ 



1 + 2{t-1)p 

□ 



Lemma 4.3 Assume now that all negative particles are free but positive particles form infinite 
life-time clusters. Then V{p) = ^ — 1). 

Proof. According to our assumptions after at most r + 1 time steps the leading particles in 
the infinite life-time cluster will take part only in short interactions with the root of the cluster. 
Consider now the flux of positive particles through the root of the cluster. According to the 
argument above after at most r + 1 iterations the flux become equal to 1/r (i.e. after each r 
moments of time a leading particle from the cluster body exchanges positions with the cluster 
root). Thus 1/r = p{V{p, p) + V{p, p)) and hence V{p) = i (i - 1). □ 

Observe that the difl'erence between short and long interactions does not matter both in the 
free and jammed phases because in the former it does not change the timing, while in the latter 
in the steady state it cannot happen if the length of the jam is greater than 1 (there are no gaps 
between particles in the cluster body). 

The derivation of the average velocity in the free phase is based on the assumption that 
V{p, p) describes the flow of noninteracting particles. If the density of negative particles is 
becoming large enough then they start interact between themselves. Assume that this density 
is so large that the corresponding velocity belongs to the jammed region and hence 

V{p,p) = l-{T-l){V{p,p) + V{p,~p))p 

y{prp) =- (^-1). 

r p 

Then V{p,p) = 1/r which does not depend on the particle density and corresponds to the case 
when all the time a positive particle interacts with negative ones, i.e. the negative particles form 
an inflnite life-time cluster. Clearly this means that if a configuration belongs to Free!^ fl Jam^^ 
then the average velocity of positive particles is equal to 1/r. 



5 Proof of Theorems 11.21 and 11.11 

These results will be proven more or less simultaneously through a series of technical Lemmas. 



Lemma 5.1 Let a configuration x satisfies the condition 



{T + 2)p{x)-{T-l)~p{x)<l. (5.1) 

Then only finite life-time (proto) clusters may he present in the configuration x. Moreover, 
there exists a partition of the integer lattice into nonoverlapping BAs of these finite life-time 
(proto) clusters and their complements. 

Proof. Assume on the contrary that there is a (proto)cluster J = x[m',m] with an infinite 
BA. By Theorem 13.11 and the definition of the functional W for any integer n > 1 we have 
VF(a;[m — n,m]) > 0. Rewriting this inequality in terms of particle densities we get 

<l{xm) +(r + l)np(x, [m — n,m — 1]) — Tnp{x, [m — n,m — 1]) 
— n(l — p{x, [m — n,m — 1]) — p{x, [m — n,m — 1])) > 0. 

Dividing by n and passing to the hmit as n — > oo we obtain 

{T + 2)p{x)-{T-l)p{x)>l, 

which contradicts to the assumption (jS.ip . 

To prove the second claim we need to demonstrate that the finite BAs cannot be infinitely 
nested. The latter means that there exists a sequence of (proto) clusters {Jn}n and an integer 
N such that for any n > N we have BA(J„) DBA(JAr). If this would be the case a positive 
particle located initially to the right of Jn will never become free despite the amount of time it 
spends in each of the clusters is finite. 

Assume that such an infinitely nested sequence of (proto) clusters {Jn}n does exist. Then 
for some m E Z and any integer n > 1 by Theorem 13. II we have W{x[m,m -\- n]) > 0. Rewriting 
this in terms of particle densities we get 

q{xm+n) +{t + l)np{x, [m, m + n — 1]) — Tnp{x, [m, m + n — 1]) 
— n(l — p{x, [m, m + n — 1]) — p{x, [m, m + n — 1])) > 0. 

Exactly as in the previous case we divide this inequality by n and after passing to the limit as 
n — >■ oo we come to the contradiction to the assumption (jS.ip . □ 

Remark. The importance of the 2nd claim of Lemma 15.11 is that it might be possible that all 
clusters have only finite life-times but a particle is joining them one after another and is never 
becoming free. 

Lemma 5.2 Let a configuration x satisfy the conditions of Lemma 15. il Then any positive 
particle after some initial transient period will stop interacting with other particles of the same 
sign (i.e. will become free). 

Proof. By the previous Lemma there exists a partition of the integer lattice into nonoverlapping 
BAs of finite life-time (proto)clusters and their complements. Choose one of these BAs and 
consider a particle ^ located originally immediately before this BA. By the definition of the BA 
the particle ^ does not belong to any BA of (proto) clusters located ahead of it and hence it will 
never join any cluster of positive particles. □ 

Results of Lemmas 15.11 15.21 prove claims (a,b) of Theorem 11.11 



Lemma 5.3 Let a configuration x satisfy the condition 

{r + l)p{x)-{T-l)~p{x)>l. (5.2) 

Then for each N ^ "L there is an infinite life-time (proto) cluster in x located to the right from 
the site N . 

Proof. Fix some iV G Z and assume on the contrary that any positive particle in x located at 
a site i > N is either free or belongs to a finite life-time (proto) cluster. 

Thus there exists a sequence of finite BAs of (proto) clusters such that BA{Jk) 5 [jk,n^k] 
and j'j^ < ruk < j'^j^i < mk+i '^^^^ oo. By the definition of for any £ > 1 we have 

e 

Y,W'{x[j',,mk])<0. 

k=l 

Between the segments there might be some additional positive particles not belonging 

to any BA and thus each of them satisfies the condition that either the number of vacancies 
immediately ahead of it exceeds r or it is preceded by a negative particle. In both cases (as 
well if there are additional negative particles) the value of the functional W' on the complete 
segment [j[,mi] instead of the union of segments is again negative. 

Denoting m := j[, := — j[ and rewriting the inequality W'{x[m,m + — 1]) < in 
terms of particle densities we get 

q{xm+ne) +(t + l)np{x, [m, m + Ui - 1]) - Tnp{x, [m, m + n^- 1]) 
— n(l — p{x^ [m, m + n£ — 1]) — p{x, [m, m + n^ — 1])) < 0. 

Dividing this inequality by n and passing to the limit as ^ ^ oo we obtain 

(r + l)p(x)-(r-l)p(x)<l. 

Thus we came to the contradiction to the assumption ()5.2p . □ 

Claims (c,d) of Theorem 11.11 follow from this result. 

In the hysteresis phase the protoclusters might be present, which makes the main difference 
to the free phase. Moreover, it might be possible that despite the protoclusters have only 
finite life-times their BAs are infinitely nested which prevents positive particles to become free 
eventually. 

Now we are ready to prove Theorem 11.21 

We start with the case {p{x), p{x)) € A. Then the particle densities satisfy the system of 
inequalities 

(r + 2)p< l + (r-l)p, 
{t + 2)p< 1 + (t-1)/>. 

By Lemma 15.21 and the first of these inequalities all positive particles will eventually become 
free. On the other hand, applying the duality relation together with the second inequality, by 
Lemma 15.21 we get the same property for negative particles. Therefore exact relations for the 
average velocities {V{x),V{x)) follow from Lemma 14.21 
Assume now that {p{x), p{x)) G C. Then 



(t + 1)p>1 + {t-1)p, 



(r + 2)p< l + (T-l)p. 



The first of these inequahties implies by Lemma [5.3l that in the configuration x there are infinitely 
many infinite life-time (proto)clusters of positive particles, while from the second inequality it 
follows that all negative particles will eventually become free. Thus some negative particle will 
become a short root of an infinite life-time cluster and from that moment of time will move 
exactly by one site in r time steps (exchanging positions with positive particles from the body 
of the cluster). Now since the average velocity is well defined for a single particle, by Lemma [2.4l 
the same result holds for all particles of the same sign and thus V{x) = 1/r. Once we get this 
we obtain the relation for the average velocity for positive particles by the duality relation. 

The case {p{x), p{x)) £ B is considered similarly to the previous one, except for in this case 
positive particles get the average velocity 1/r. 

Theorems 11.21 and 11.11 are proven. □ 



6 Hysteresis region 

Our previous results are concerned with the behavior of configurations in "pure" free and jammed 
phases. Now we shall show that in the region H located between these phases (see FigHJ^right)) 
the coexistence of free and jammed configurations takes place. In other words in this region 
there might be both freely moving and jammed configurations with the same particle densities. 
Due to the duality between positive and negative particles in what follows we consider only the 
case when p > p. 

Lemma 6.1 Configurations from the set PreeJ^ are dense in := {{p,p) £ H : {t + l)p < 
l + (r-l)p}. 

Proof. Consider a family of spatially periodic configurations parameterized by nonnegative 
integers n,m, ki, . . . , kn (only one spatial period is shown): 

(r+l)n+2iV 2 m 



+ ^ + ^...+ +^...+^ (6.1) 

T+2ki T+2fc2 r+2fc„ 

Setting := ki the length of the period can be written as (r -|- l)n -|- 2(m -|- N). The value 

r -|- 2ki is equal to the duration of the short interaction plus an arbitrary even number. This 
guarantees that only short interactions may take place under dynamics. On the other hand, the 
number of vacancies between consecutive positive particles is not smaller than r, which implies 
that for any ki a positive particle cannot catch up with another positive particle (since only 
short interactions may take place). Therefore x £ Free+. 
Let a configuration x belongs to the family (j6.ip . Then 

r \ n + m ( at\ 

P^^^ = (r + l)n + 2(m + iV) =^"^"'"^'^^' 

, m _ 

p(x) = -, r ; — =: r{n,m, I\ ). 

' {t + l)n + 2{m + N) ^ ' ' ^ 

Observe now that for any u,v^z >{) such that u + v > Q the identities 

r{u,v,z) =r( — — , — — , — — , r{u,v,z) =r — — , — — , — — 
u + v u + v u + v u + v u + v u + v 



hold true and that 



Therefore the "right" boundary of the region H (corresponding to the case z = 0) contains 
configurations from Free+. Moreover, the points on the "right" boundary for which the above 
representation takes place are dense. Additionally we see that these functions decay with respect 
to each their argument. Therefore for a given pair or real numbers u,v > 0, u + v > the 
pair of functions (r{u,v,z), r{u,v,z)), considered as functions of the third real argument z, 
parameterizes a curve Ru^v{z) '■= {r{u,v,z), f{u,v,z)) which starts (when z = 0) at some point 
on the "right" boundary of H and goes through H until it reaches the "left" boundary. 
For d > 1 denote £>d((ni, . . .,Ud), iu[, . . .,u'^)) := Y^i " K\- Then 

( \ -D I i\\ ^ 2rg3((^,^;,2:),(^^^^zO) 

Q2KRuA^),Ru',v'{.Z ) < - — ■ ■ — 6.2) 

(u + f + z) [u +v' + z') 

for any pair of triples (u, f , z), («', u', z'). Therefore the curves Ru,v with u,v £ M+ U {0} fill in 
the entire region H-^.. 

It remains to show that for any real (p, p) € and any e > there exists a triple of 
nonnegative integers {n,m, N) such that 

Q2{{p,p), RnM^)) < 

To do this set n' := m' := Then for any K > the curve Rn' ,m' goes through the point 
(p, p) at some real parameter value N' . Since points on the "right" boundary of H represented 
by integer values of (n, m, A^) are dense we may choose two sequences of nonnegative integers 
{rifc}, {mk} ^^^^ oo such that 



Therefore setting 



Q2{Rn',m'{0), i?„fe,mfe(0)) ''-^ 0. 



p , Kp , N' 
u := -, V := -, z — 



i~' i~' /i /' 

p + p P + P n' + m' 

rik mi. N 

Uk := ■ , Vk := ■ , Zk :-- 



we get 

Using ()6.2p we obtain 



nk + mk nk + ruk Uk + 



Q2[{uk,Vk),{u\v')) 0. 



( \ u ( '\\ ^ '^Tg3{{uk,Vk,Zk),{u',v',z' 



(1 + Zfc)(l + Z') 
< 2T{Q2{{Uk,Vk), {u',v')) + \Zk - z'\). 



Therefore the curves Rnk,mk come arbitrary close to the point {p,p) as A; — > oo at some integer 
parameter values N, for which 

I / 1 k—*oo „ 

\zk- z\ — > 

and this finishes the proof. □ 



To study the "left" boundary of H consider another family of spatially periodic configura- 
tions: 

(T+2)n+2N 3m 



+ ^ + ^ •••+ ^ i^i^-'-i^' 

T+l+2fci T+l+2fc2 T+l+2k„ 



The notation here is the same as for the family (|6.ip and the length of the period is equal 
to (r + 2)n + 2(m + A^). A direct check shows that again in this family both positive and 
negative particles are free. The difference is the absence of long interactions under dynamics 
for configurations from (j6.ip and their inevitable presence for configurations from ()6.3p for an 
arbitrary large time. If a configuration x belongs to the family (j6.3p then 

n + m _ m 

Py^) = 7 — rrr^ — tttt — Py^) 



{t + 2)n + 2{m + N) ' (r + 2)n + 2(m + AT) 

Setting again = we obtain the "left" boundary of ff, namely 

(r + 2)p(x)^l + (T-l)p(x). 
On the other hand, for the first family 

l>(p-p)(T + l) + 2p =^ p<i±^(Lpl) 

r + 1 

while for the second family 

l + /5(r-l) 



1> (p-p)(t + 2) + 3/5 =^ p< 



r + 2 



Comparing these inequalities with the critical densities p'^ and pc we see that the family ()6.3 
belongs to the free phase while the family (16. ip does not belong there. 



It remains to show that in this region there are jammed (for all times) configurations as well. 
According to Lemma IS.ll the inequality (j5.ip (defining the "right" boundary of Hj^) guarantees 
the absence of infinite life-time (proto) clusters in H. This, together with the denseness of 
densities corresponding to free configurations in H (see Lemma 16. ip explains the subtlety of 
this problem: one needs to find configurations having only (proto) clusters of arbitrary large but 
finite life-time. To understand how this can happen consider the dynamics of a spatially periodic 
configuration x G X with r = 2: 

lOOOlili 
01002222 
OOlOilil 
10202220 
01303110 

olioiiol 

12202200 

liioiloo 

22202010 
11303001 
21101002 
12200101 



Using the notation introduced in the proof of Lemma 16.11 we have n = 3,m = 2,A^ = 3. 
Thus r{n,m, N) = 3/8, r{n,m, N) = 1/4 and hence 

(r + l)r(n, m, A^) = 9/8 < 1 + (r - l)f(n, m, N) = 5/4 < (r + 2)r(n, m, N) = 3/2, 

i.e. X belongs to the region H. During the dynamics each positive particle in makes short 

and long interactions with negative ones. Finite (proto) clusters are present for all time. On the 
other hand only trivial life-time clusters are present (at moments of time when one of positive 
particles catches up with another), hence the configuration in this example does not belong to 
the jammed phase. Nevertheless at any moment of time there are infinitely nested BAs of trivial 
protoclusters. 

Our numerical simulations, which were carried out for a broad region of spatial periods L 
(up to L = 200) indicate that for each spatial period L for all pairs of positive integers n, m for 
which n + m < L and (r(n, m, L — n — m),r{n, m, L — n — m)) G H there were configurations 
with only trivial (proto)clusters for arbitrary large time. This allows to formulate the following 

Hypothesis. For each pair of densities belonging to the region H there exists an admissible 
configuration having positive and negative particles with these densities for which only trivial 
(proto) clusters are present for arbitrary large time. 

For spatially periodic configurations the existence of average velocities follows from the obvi- 
ous time-periodicity of the dynamics of particles. However in the general case this is not proven. 
Moreover, without the assumption about the existence of particle densities one can construct 
configurations for which average velocities do not exist. 

7 Dynamics of an active tracer 

To model the motion of an active tracer in the traffic flow we add to the flow a single particle 
(representing the tracer). If the tracer moves in the same direction as the flow we think about 
this particle as having the same sign as particles of the flow, and having the opposite sign when 
the tracer moves against the fiow. In distinction to the passive tracer studied in [5j the active 
one changes the behavior of the flow due to interactions with other particles. Clearly a single 
positive particle added to the fiow of positive particles do not change its asymptotic behavior. 
On the other hand the case when the tracer moves against the flow is equivalent to the model 
discussed in this paper with zero density of positive particles (tracer) and a certain density of 
negative particles (the flow) . According to Theorem 11.21 we get the following average velocity of 
the active tracer as a function of the density of the fiow p: 



Here the average velocity V+(/3) corresponds to the case when the tracer moves in the same 
direction as the fiow, while V-{p) describes the case when the tracer moves in the opposite 
direction. Fig. [2] shows that there is a region of fiow densities when it is more advantageous to 
move against the fiow than along it. The explanation is that in high density region the time 
necessary to exchange positions with a particle moving in the opposite direction is becoming 
smaller than the time a particle from the flow spends in inevitably being present large clusters. 

^The sequence of types of interactions for a single particle is periodic but may be arbitrary. 
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Figure 2: Dependence of the average velocity of the active tracer V± on the particle density p. 
V+ - along the flow (thick line), V- - against it (thin line). 

8 Configurations having no particle densities 

So far we were discussing statistical properties of configurations for which both positive and 
negative particle densities are well defined. In general one or both of these densities might be 
not well defined and one needs to deal with lower and upper densities. Recall that Lemma |2. II 
claims their invariance with respect to dynamics. The following result extends the description 
of the Phase diagram for this more general case. 

Theorem 8.1 Let x G X 

(a) {t + 2)p+{x) < 1 + (r 

(b) (r + 2)p+(x) < 1 + (t 

(c) {T + l)p-{x) > 1 + (r 

(d) {t + 1)p~{x) > l + (r 

The proof is similar to the proof of Theorem 1 1 . 1 1 and therefore we skip it. 

In distinction to the case when the densities are well defined, here we cannot prove the 
existence of average velocities. Instead one considers upper and lower particle velocities for 
which it is possible to get a representation of the same sort as in Theorem 11.21 
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